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Stripe width variation in ultrathin magnetic films is a well known phenomenon still not well 
understood. We analyze this problem considering a 2D Heisenberg model with ferromagnetic ex- 
change interactions, dipolar interactions and perpendicular anisotropy, relevant e.g. in Fe/Cu(001) 
films. By extending a classic result of Yafet & Gyorgy (YG) and using Monte Carlo simulations 
we calculate the complete zero temperature phase diagram of the model. Through this calculation 
we analyze the correlation between domain walls structure and stripe width variation, as the per- 
pendicular anisotropy changes. In particular, we found evidences that the recently detected canted 
state becomes the ground state of the system close to the Spin Reorientation Transition (SRT) 
for any value of the exchange to dipolar couplings ratio. Far away of the SRT the canted ground 
state is replaced by a saturated stripes state, in which in-plane magnetization components are only 
present inside the walls. We find that the domain wall structure strongly depends on the perpen- 
dicular anisotropy: close to SRT it is well described by YG approximation, but a strong departure 
is observed in the large anisotropy limit. Moreover, we show that stripe width variation is directly 
related to domain wall width variation with the anisotropy. 

PACS numbers: 75.40. Gb, 75.40.Mg, 75.10.H 
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I. INTRODUCTION 

Pattern formation in ferromagnetic thin films with 

perpendicular anisotropy and its thermodynamical 

description have been the subject of intense experimen- 

taiLaa ifiALBA ll L U ^ theoretical 12 - 1 ^ 14 - 1 ^ 16 -^. 1 ^ 19 .^ 

and numerical 2 ^^^^^^' 2 ^ WO rk in the last 

20 years. Magnetic order in ultrathin ferromagnetic 

films is very complex due to the competition between 

exchange and dipolar interactions on different length 

scales, together with a strong influence of shape and 

magnetocrystalline anisotropies of the sample. These in 

turn are very susceptible to the growth conditions of the 
films 8 i 3i_ 

Among the different magnetization patterns that have 
been observed in these systems, striped order (i.e., modu- 
lated patterns of local perpendicular magnetization with 
a well defined half-wavelength or stripe width h) at low 
temperatures is an ubiquitous phenomenon. One in- 
triguing fact is the strong variation displayed by the 
equilibrium stripe width h in many of these systems, 



when either the temperature or the film thickness is 
changed^^^^. The origin of such variation is still 
controversial, but recent results suggest that a key point 
to understand it is the role played by the interfaces (i.e., 
the domain walls) between stripes^. Thus, a starting 
point to study this problem is to compare the energies 
of striped patterns with different domain wall configu- 
rations. An accurate description of the domain walls 
requires to take into account not only the perpendicu- 
lar component of the local magnetization, but also the 
in-plane component. Indeed, some experimental results^ 
are consistent with the presence of Bloch domain walls, 
as expected for perpendicularly oriented magnetization 
domains^ 2 -. 

To compute the energy contribution of domain walls 
it is important to consider explicitly the out of plane 
anisotropy, together with the exchange and dipolar in- 
teractions, whose competition is the responsible for the 
appearance of striped patterns. A minimum model that 
contains all these ingredients is the 2D dimensionless 
Heisenberg Hamiltonian: 
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where Si are classical unit vectors, the exchange and 
anisotropy constants are normalized relative to the dipo- 
lar coupling constant (6 = J/0, 77 = K/Sl), < i,j > 
stands for a sum over nearest neighbors pairs of sites in 
a square lattice, stands for a sum over all distinct 

pairs and nj = \n — fj\ is the distance between spins i 
and j. 

In the large anisotropy limit 77 — ► oo this model reduces 
to an Ising model with short range ferromagnetic and 
long range antiferromagnetic interactions, whose ground 
state is the striped one^. In that limit the stripe 
width increases exponentially with the exchange to dipo- 
lar coupling rati o 21 ' 26 h ~ exp(<5/2). For low values of 
the anisotropy, the ground state of this model changes 
to a planar ferromagnetic stated. In a classic work, 
Yafet and Gyorgy computed the energy of striped do- 
main configurations with Bloch domain walls, by con- 
sidering a sinusoidal (perpendicular) magnetization pro- 
file at the walls and saturated magnetization inside the 
stripes^. They found that above certain threshold value 
V > f]min a striped configuration has less energy than 
a uniformly (in-plane) magnetized one. At this point 
the system shows a Spin Reorientation Transition (SRT) . 
Also within this approximation the stripe width shows an 
exponential increase with the anisotropy strength, while 
the domain wall width decreases algebraically for large 
enough values of 77, at least for large values of 5, This 
approximation is expected to work well close to the SRT, 
where the effective anisotropy is smalU&. However, for 
large values of 5 and 77 not too close to the SRT (i.e., 
when the width of both stripes and walls are large) one 
can expect a large departure in the wall energy contribu- 
tion with respect to the true domain wall configuration, 
which should approach a hyperbolic tangent profile^ 2 -. 

Another important point concerns the in-plane mag- 
netization component inside the stripes. In their work 
Yafet and Gyorgy considered striped solutions where the 
only in-plane components lay inside the walls, although 
they pointed out how to extend their calculations to con- 
sider non-saturated magnetization inside the stripes^. 
Based on that calculation, Politic reported that, at least 
for large enough values of 5, the magnetization should 
show an abrupt saturation very close to the SRT, sug- 
gesting that the in-plane component inside the stripes is 
not relevant. However, recently Whitehead et ob- 
tained evidences of a non-saturated ground state at a 
relatively small value of S for a wide range of values of 
the anisotropy strength 77. Their numerical results also 
showed a correlation between the stripe width and the in- 
plane component variations and suggest that this ground 
state configuration stabilizes at finite temperature, giv- 
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ing rise to what they called a "canted" phase. Hence, it is 
important to revise the zero temperature phase diagram 
of this model in the whole (S, 77) space, including canted 
configurations, in order to determine to what extent they 
can be relevant to real systems and their possible influ- 
ence to the stripe width variation phenomenon. 

It is worth noting that the perpendicular anisotropy 
changes (inversely) with the film thickness 2 ^.. It has also 
been pointed out that the changes in the film thick- 
ness can act as a change in the effective temperature^. 
Hence, understanding the variation of the equilibrium 
stripe width and the associated domain wall structure 
as the anisotropy changes at zero temperature can be of 
great help to understand the corresponding properties a 
finite temperature, a far more complex problem. 

In this work we analyze the complete equilibrium phase 
diagram at zero temperature in the (8, 77) space of Hamil- 
tonian (!]). Following Yafet and Gyorgy 's work, we only 
consider straight domains, i.e., domains in which the spin 
orientation can be modulated along the x direction but is 
constant in the perpendicular direction y ((x,y) are the 
coordinates on the plane of the film ) . We also consider 
only Bloch walls, i.e., walls in which the magnetization 
stays inside the yz plane. In order to clarify notation 
and units, and to include some extensions, we briefly re- 
view Yafet & Gyorgy's approximation in Appendix[Al To 
obtain the phase diagram we compute the energy of dif- 
ferent types of magnetization profiles and compare them 
with simulation results obtained through a zero temper- 
ature Monte Carlo (MC) method specially designed for 
the present purposes. The MC method is presented in 
Appendix El In section Ull we derive a general expression 
for the energy of striped magnetization profiles. The zero 
temperature phase diagram and associated properties are 
derived in section IIIII A discussion and conclusions are 
presented in section llVl 



II. ENERGY OF STRIPED MAGNETIZATION 
PROFILES 

Let us consider a square lattice with N — L x L 
sites, characterized by the integer indexes (x, y), where 
-L/2 < x < L/2 and -L/2 < y < L/2, in the limit 
L — > 00. Hence, the index i in Eq.(U) denotes a pair 
of coordinates (x,y). Following YG approximation^, we 
consider only uniformly magnetized solutions along ev- 
ery vertical line of sites, i.e. S( XjV ) = M(x), Vy and allow 
only Bloch walls between domains of perpendicular mag- 
netization, i.e. M x (x) — Vx. Therefore, 
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YG showed that for these types of spin configurations 
the energy per spin can be mapped onto the energy of 
(2) a one dimensional XY model. The energy difference be- 
tween an arbitrary magnetization profile M[x) and a uni- 
Then, for every value of a; there is only one independent formly in-plane magnetized state is then given by (see 



M(x) = [lf(i)f + [M y {x)f = 1. 



component of the magnetization. 



Appendix IA 1 



M(x)] ={6- 2c 2 ) - - E M{x).M(x + 1) 



\x-x'\ 

x.x' 1 1 



where 6' = 8 - 2 c x , k' = r\ - 3g, ci = 0.01243 .. ., c 2 = 
0.07276 g= 1.202057 ... and 

C = C[M v (x)] = 2(c 2 -ci)-^2M y (x)M v (x+l) (4) 



Although small, this correction term makes a non negli- 
gible contribution when the domain walls are of the same 
order of the lattice constant. This happens for small val- 
ues of S (5 < 5), where both the stripe and wall widths 
are of the order of a few lattice spacings. For larger values 
of S it is reasonable to assume a smooth magnetization 
profile^ M y (x + 1) « M y (x), so that 



C m 2(c 2 - ci 



1 - (M z (x)Y 



(5) 



can be absorbed into the anisotropy term in Eq.j3]), re- 
placing k'-^k — rj — 3g + 2(c 2 — ci). 

Considering now a stripe-like periodic structure of the 
magnetization profile with period 2h, i.e. 



M z (x + h) = -M z (x), 
we can make use of a Fourier expansion: 

M z (x) = M E b ™ cos 



(6) 



(7) 



where we have assumed M z (x) an even function of x just 
for simplicity. 

The anisotropy term in Eq.j3]) can be written as 



, 2h 



and the dipolar term^ 
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In the general case, it is better to let the exchange term 
(and the correction term C as well) in Eq. |(3j) expressed 
in terms of the angle 4>(x) between M{x) and the z axis: 



M z (x) = cos[(j)(x)} 
M v (x) = sm[<f>(x)] 



(11) 
(12) 



where the angle (j)(x) has the same periodicity of M(x). 
We have that 



= cos Ma)- 4>{x+l)}. (13) 



Putting all the terms together we get the general ex- 
pression: 



M;8, V ] = (5-2c 2 )-5'-Ecos[ 



[X) - 0{X 



k'M? 



1)]+M 2 E b 2 m D m (h)-^ E b2 m + C - ( 14 ) 



m— 1,3,... 



m — 1,3,. 



III. ZERO TEMPERATURE PHASE DIAGRAM magnetization profiles M z {x) and compare the energies 



In this section we look for the minimum of Eq. (fl4|) 
for different values of S, r\. We propose different striped 
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obtained by minimizing Eq. (fT4|) for each profile with re- 
spect to variational parameters. 

A. Small values of 8: Sinusoidal Wall 
magnetization Profile (SWP) approximation 

We first consider a profile as proposed by YG, that is 
constant inside the stripes and has a sinusoidal variation 



where Mq is the absolute value of magnetization inside 
the stripes and w is the wall width. In order to allow 
for canted profiles, we take Mo — cos 0, where 9 is the 
canted angle, i.e. we define it as the minimum angle of 
the local magnetization with respect to the z axis. Yafet 
and Gyorgy solved this variational problem for Mq = 1 
in the continuum limit^, i.e. when h ^> 1 and w 1, 
so that the profile can be considered a smooth function 
of x. While this approximation is expected to work well 
for large enough values of 5, it breaks down for relatively 
small values of it, where the discrete character of the 
lattice has to be taken into account. However, the varia- 
tional problem for that range of values of S can be solved 
exactly (although numerically) by minimizing Eq. (fT4|) 
with respect to the integer variational parameters h and 
w and continuous parameter 9. In other words, for every 
pair of values (6, if) we evaluate the energy Eq. (fl4|) for the 
profile lfl~5|) with different combinations of h — 1,2,... 
and w = 1,2,... within a limited set. For every pair 
of values h, w,we look for the value of 9 that minimizes 
the energy with a resolution A9 = 0.01 and compare all 
those energies. Details of that evaluation are given in 
Appendix [B] This calculation is feasible for values up to 
S = 10, for which the maximum value of h (bounded by 
the stripe width in the 77-^00 limit) remains relatively 
small (smaller than h = 140). Some results for S = 12 
close to the SRT were also obtained. All the results of 
this calculation are compared against Monte Carlo (MC) 
simulations. Details of the MC method used are given 
in Appendix El Through these calculations we obtain a 
zero temperature phase diagram for low values of S. 

Before presenting the results, it is important to intro- 
duce some notations and definitions of different types of 
solutions. We distinguish between four types of solutions. 
If the minimum energy solution corresponds to w = 1 and 
9 = (within the resolution AO), we call this a Striped 
Ising Profile (SIP), i.e. a square wave like profile. If 9 = 
but w > 1, we call this a Saturated State. These states 
show a finite parallel component of the magnetization in- 



inside the walls between stripes (see Fig.l in Ref. [13): 



(15) 
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side the walls. IfO < 9 < ir/2 the solution is a Canted 
State. Finally, if 6* = ir/2 (Mq = 0) we have a Planar 
Ferromagnet (PF). 

The zero temperature phase diagram for small values of 
5 (6 < 5) is shown in FigEJ For relatively large values of 
i] the minimum energy configuration is always the Ising 
one (SIP), with a stripe width independent of rj. For 
small values of r\ the minimum energy configuration is 
the PF, with a spin reorientation transition line (SRT), 
either to the Ising state for h < 3 (S ~ 2)) or to a canted 
one for h > 3 (S > 2). No Saturated State configurations 
are observed for 5 < 6. 

Inside the canted region, a strong stripe width varia- 
tion with the anisotropy is observed at constant 8 . Note 
that the vertical lines that separate Ising striped states 
with consecutive values of h bend inside the canted re- 
gion and become almost horizontal as 8 increases. Hence, 
the exponential increase of h with 8 in the Ising region 
(vertical lines) changes to an exponential increase with r] 
deep inside the canted region (curved lines on the right 
ofFigd). 

We also find an excellent agreement between the si- 
nusoidal wall profile approximation and the MC results, 
except close to the transition between the Ising and the 
canted states. Such disagreement is due to the fact that 
the actual wall is not well described by a sinusoidal pro- 
file far away of the SRT line, as will be shown later. In 
Fig|2] we show a comparison between the energy of the 
SWP and the MC results as a function of 77 for 8 = 4.58. 
The range of values of 77 where the walls are not well 
described by a sinusoidal profile increases with 5. 

For large enough values of S the variational problem for 
the SWP can be solved in a continuum approximation 
introduced by YG^. This leads to the equations (see 
Appendix I A 2ft : 
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FIG. 1: (Color online) Zero temperature phase diagram for 
small values of 5. Black filled symbols and black solid lines: 
MC simulations. Open red symbols: SWP approximation. 
Squares and continuous black lines correspond to transition 
lines between striped states of different width. The shaded 
region corresponds to the Canted State (0 < 9 < 7r/2). 
Triangles are transition lines between Planar Ferromagnet 
(9 — 7r/2) and Canted States (Spin Reorientation Transition 
line). Circles mark transitions between the Canted and the 
Stripes Ising state (9 = and w = 1). Notice the excellent 
agreement between the MC and SWP calculations close to the 
SRT, while the SWP approach underestimates the transition 
line between the canted and Ising Stripes states. The dotted 
line corresponds to the contimuum approximation of YG for 
the SRT (Eq.(H9|)). 
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where A = w/h, k = n/h and 7 = 7r 2 /3 — k. In the limit 
A — * 1 (pure sinusiodal profile) these equations can be 
solved analytically and predict a SRT at the line 



Vsrt{8) = a - — 



(19) 



with a = tt 2 /3 + 3.g- 2 (c 2 - a) (see Appendix |A~2|) . The 
line Eq. l|19p is also depicted in Fig[TJ Notice the disagree- 
ment between the continuum approximation and the ex- 
act one for S < 5. This discrepancy becomes smaller than 
1% only for 6 > 7. 

For arbitrary values of r\ and 6 Eqs. (fl6|) - (fl8)) can be 
solved numerically. In Fig|3] we show the numerical so- 
lutions for 9 and h as a function of 77 for different values 
of S. We see that the range of values of the anisotropy 77 
for which the canted angle is appreciable different from 
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FIG. 2: (Color online) Energy per spin (with respect to the 
PF state) vs. 7? for 8 = 4.58 within the SWP (full red line) 
and MC calculations (filled black circles). Open black circles 
correspond to the energy of a Striped Ising Profile (SIP) with 
h = 9 (equilibrium value for r\ — > 00). 



zero within this approximation is strongly depressed as 
8 increases. For values 6 ~ 100 the canted configuration 
almost disappears, except very close to the reorientation 
line, consistently with the results reported by Politic. 

Indeed, from our MC simulations, we observe that the 
range of values for which the canted state has the min- 
imum energy gradually shrinks as 6 increases, being re- 
placed by a saturated state for values of r\ above certain 
threshold. This can be observed in FigJH where we show 
the behavior of the canted angle and the in-plane magne- 
tization component M\\ = (1/L) J2 X M v (x) as a function 
of 77 for S = 7.5. The Monte Carlo data shows the exis- 
tence of a wide range of values of 77 for which the canted 
angle is zero while Mm ^ 0, meaning that the non null in- 
plane components are concentrated inside the walls. In 
other words, in that region we have a saturated state with 
thick walls w > 1. Notice also that the SWP approach 
completely fails to describe those states. Moreover, we 
observe from our MC simulations that the SWP cease 
to be the minimum energy solution for values of 77 rela- 
tively close to the SRT, well before the saturated state 
sets up (see Figj4j) . This effect becomes more marked as 
S increases. 



B. Large values of S: Hyperbolic Wall 
magnetization Profile (HWP) approximation 

As already pointed out, the actual magnetization pro- 
file departs from the SWP for large values of 77 and 5. 
This is expected from micromagnetic theory, which in 
that limit predicts that the wall structure will be domi- 
nated by the interplay between anisotropy and exchange, 
leading to an hyperbolic tangent shape of the wall^. This 
can be observed in FigUl Hence we considered a peri- 
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FIG. 3: (Color online) (a) Canted angle 8 and (b) equilibrium 
stripe widht h, as a function of the anisotropy r\ for different 
values of S from the continuum approximation of the SWP 
Eqs. l|16ll - I|18p . The full grey line at the left corresponds to the 
SRT borderline given by Eq. lfl9ll with h = S (see Appendix 
lA2i 



the large 5 limit, assuming a smooth profile h 3> 1 and 
l w 1, the anisotropy energy can be expressed as: 
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FIG. 4: (Color online) Canted angle (circles) and in-plane 
magnetization (squares) as a function of 77 for S — 7.5. Filled 
black symbols correspond to MC calculation. Open red sym- 
bols corresponds to the discrete SWP approximation, while 
the red doted and dashed lines correspond to the continuum 
(YG) approximation of the SWP. Continuous black and red 
lines are only a guide to the eye. 



odic magnetization profile with hyperbolic tangent walls 
(HWP) defined, for a wall centered at x — 0, by 



M z (x) = M tanh 



for -h/2 < x < h/2, (20) 



hl w Jo 



d<p{x) 
dx 
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Mn 



tanh u 



du. (22) 



together with Eq.|(6]) where Mq — cosO as before. In Solving the last integral we finally obtain 



= -5 



h \ M 



Mo tanh 
I 



h 

w 

— 111 



tanh ( — — 



(23) 



The dipolar energy can be calculated using Eqs.© and 
(flOll . The Fourier coefficients for the profile ff20|) can be 
computed using the approximation (|D1|) (see Appendix 
ID]) . This leads to an expression for the total energy as a 



function of the variational parameters h, 8 and l w that 
can be minimized numerically. Comparing the minimum 
energy solution for the SWP and the HWP we obtain 
the crossover line between sinusoidal and hyperbolic wall 
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structure shown in FiglB] (dashed line). Above that line 
the HWP has always less energy than the SWP. We also 
calculated the transition line between the canted and the 
saturated states by setting the condition 6 = 0.01, to 
be consistent with the criterium used in the MC calcula- 
tions. The results are shown in Fig|6] together with the 
SRT line Eq. (fl9|) . and compared with MC calculations 
up to 5 — 15. The excellent agreement with the MC 
results gives support to the analytic approximations. 

For large values of r\ the exponential increase of h 
makes it cumbersome to apply the approximation of Ap- 
pendix [D] for the calculation of the dipolar energy. In- 
stead of that, we can use the following heuristic argument 
to obtain a reasonable approximation. The main error in- 
troduced by the SWP approach is in the exchange and 
anisotropy contributions to the energy. Since the main 
contribution to the dipolar energy is given by the inter- 
action between domains, we can assume that the dipolar 
contribution of the wall is relatively independent of its 
shape. Hence, we can approximate it by Eq. (|A13| l. Fur- 
thermore, taking w — f l w (/ is a fitting parameter of or- 
der one to be fixed later) in the limit A < 1 C 1), 
G(A) is very well approximated by^ 



M\x) 
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(25) 



we compare the energy obtained with the above equation 
with that obtained using the approximation of Appendix 
[D]for different values of the system parameters. We ver- 
ified that the error made by the approximation Eq. ([25|) 
taking / = 4 is always smaller than 1% for h/l w > 20. 
We also observe that the best agreement with the MC re- 
sults is obtained for / = 4. Assuming then M — 1, the 
total energy per spin (relative to the parallel magnetized 
state) for the HWP can then be approached by 



ghwp — 



5/l w + 21 w (k-tt 2 /3) 



3 nh 



10 L 



(26) 



Minimizing Eq. ([26"l) with respect to the variational pa- 
rameters h and l w (using Eq. (|24")l ) leads to: 
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FIG. 5: (Color online) Magnetization profile in a saturated 
state (9 = n/2 and w > 1) for 5 = 8 and r\ = 7. Black symbols 
correspond to the MC simulations (M = 0.999 and h = 20). 
The full red line corresponds to the YG approximation for the 
SWP (M = 0.98, w = 6.7 and h = 20). The black dashed 
line is a fit using a hyperbolic tangent function Mo tanh(x/l m ) 
(M = 0.993 and l w = 1.48). 
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FIG. 6: (Color online) Zero temperature phase diagram for 
large values of 5. The shaded region corresponds to the 
Canted state. Symbols correspond to MC simulations and 
lines to theoretical results. The dashed line correspond to the 
crossover between sinusoidal and hyperbolic wall structure. 
The lower line (blue) corresponds to Eg. I(19[) . The middle 
line (red) is obtained from the HWP minimum energy solu- 
tion with 8 — 0.01. The upper line (black) corresponds to 
Eq,p3]l. 
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in agreement with a derivation made by Politi 



With the previous calculation we can also estimate 
the transition line between the saturated and the Ising 
Striped state. In the large h limit the energy for a SIP, 
i.e. for 
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can be easily calculated from Eg • (|14|) . The Fourier coef- 
ficients can be obtained as the A — » limit of Eq. (|All|) : 



Using Eq. lflO]) the dipolar energy is then given by 



(30) 
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where /3 = 7 e + In 4 - 1, 7 e « 0.577216 is the Euler 
gamma constant and ip(x) is the digamma function^. 
The energy per spin respect to the in-plane magnetized 
state is then given by 



7T 2 25'- (3 4ip(h) 



(32) 



Minimizing Eq. (f32|l with respect to h leads to the equa- 
tion 8'/2 - P = F(h), where F(h) = ifj{h) - hif>'(h) ~ 
In h— 1, thus recovering the known result h ~ e 5 / 2 . Com- 
paring the energies, we find that the HWP has less en- 
ergy than the Ising state for any value of 77. Eq. ([27j) 
shows that the stripe width variation in the Saturated 
state is determined by the change in the wall width as 
the anisotropy increases. Hence, h will change until the 
wall width reaches the atomic limit, i.e. for l w — 1, where 
Eq . (|27)1 recovers the Ising behavior h 



s 5 / 2 . Imposing 
the condition l w = 1 to Eq . (|28)1 we obtain the transition 
line between the Saturated and the Ising Stripes states: 

1 7T 2 

V=-6-2+ — + 3g-2(c 2 - Cl ), (33) 

which is also shown in FigEl in complete agreement with 
the MC results. 

In Fig[7]we show the variation of stripe width h vs. r\ 
for different values of 5 in the HWP, comparing the varia- 
tional solution from Egs- lpT]) and lf23j) using the approxi- 
mation of Appendix[D]and the asymptotic approximation 
given by Eqs. l(27j) and f28|) . In Fig [8] we compare the equi- 
librium stripe width h as a function of r] obtained within 
the different approximations used in this work for S = 10 
and with the MC simulations. Notice that the asymp- 
totic approximation for the HWP given by Eqs. (|27| and 
(|28)l shows a better agreement with the MC results than 
using the approximation (|D1(1 for the Fourier coefficients 
in the dipolar energy. This is because we adjusted the 
fitting parameter / to optimize the agreement with the 
MC results at low values of 6. From Fig[7] we see that 
the discrepancy between both approximations becomes 
negligible in the large 5 limit. 



500 




FIG. 7: (Color online) HWP equilibrium stripe width h vs. 
r\ for different values of 5. Full green lines correspond to the 
variational solution of Eqs. lpTjl and l[23|l using the approxi- 
mation l|Dip for the Fourier coefficients in the dipolar energy. 
The dotted lines correspond to the asymptotic approximation 
given by Eqs. (|27|) and (|28 j> . The dash-dotted line corresponds 
to the SRT borderline given by Eq-lUnj with h = S. The 
dashed line corresponds to the borderline between Saturated 
and Ising Stripes states. 
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FIG. 8: (Color online) Comparison of the equilibrium stripe 
width h vs. 77 obtained within the different methods for 
5 — 10. Symbols correspond to MC simulations. Full black 
line corresponds to the asymptotic approximation for the 
HWP given by Eqs. l(27|l and (f28|) . The red dashed line cor- 
responds to the variational solution of Eqs. lf2T|) and ll23ll us- 
ing the approximation (jD 1 j) for the Fourier coefficients in the 
dipolar energy. The blue dotted line corresponds to the con- 
tinuous SWP. The vertical dashed line corresponds to the 
transition between Saturated and Ising Stripes states given 
by Eq. (f33]> . 
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IV. DISCUSSION AND CONCLUSIONS 

The main results of this work are summarized in FigsQ] 
and [H which display the complete zero temperature 
phase diagram of the model defined by the Hamilto- 
nian JT]). Working upon reasonable assumptions for the 
ground states, like perfectly straight modulations in one 
dimension and Bloch domain walls, we analyzed mini- 
mum energy configurations combining a variational anal- 
ysis with Monte Carlo results. We find four qualita- 
tively different kinds of solutions: a planar ferromagnet 
for small anisotropies, and three types of perpendicular 
striped states: a canted state where the local magnetiza- 
tion has a finite in-plane component, a saturated state in 
which the in-plane component is restricted to the domain 
walls, and an Ising stripe state with sharp walls for large 
anisotropies. 

The canted and staturated states give valuable in- 
formation on the behaviour of the stripe width as the 
anisotropy and exchange parameters change, a still open 
and debated question^. We find that stripe width vari- 
ation is directly associated to the presence of finite width 
domain walls. For large enough values of the anisotropy rj 
the ground state of the system is always an Ising Striped 
state, no matter the value of the exchange coupling 8. In 
those states domain walls are sharp and the stripe width 
is completely independent of rj. It grows exponentially 
with the exchange coupling. 

At the SRT the system always passes through a canted 
state as the anisotropy increases, although the range 
of values of rj where the canted angle is different from 
zero narrows as 5 increases. The exchange to dipo- 
lar coupling ratio in fee Fe based ultrathin films can 
be roughly estimated to be S ~ 100 (considering a cu- 
bic bilayer of Fe/Cu(100), where 4 - the exchange coupling 
Ji? e ~ 30meV, the lattice constant c?Fe ~ 2ML and 34 
HFe ~ 3 /is)- For 5 ~ 100 the anisotropy interval for 
the canted phase is approximately A77 = rj — rjsRT ~ 0.2. 
Although narrow, this suggests that the canted phase 
should be detectable close enough of the SRT, in systems 
like low temperature grown^ Fe/Cu(100) or Fe/Ni/Cu 
films^ (room temperature grown Fe/Cu(100) do not ex- 
hibit SRT—, suggesting a rather large value of the micro- 
scopic anisotropy). 

For S < 6 stripe width variation appears always to- 
gether with a varying canted angle. Close enough to the 
SRT domain walls present a sinusoidal shape in agree- 
ment with YG results, but as the anisotropy and the ex- 
change increase, the wall profile changes to a hyperbolic 
tangent shape, as expected from micromagnetic calcu- 
lations, while the magnetization inside the domains be- 
comes fully saturated. For S > 6 the ground state is given 
by the Saturated State, except very close to the SRT. A 
similar effect (i.e. a crossover between a sinusoidal and 
a saturated magnetization profile) has been observed in 
room temperature grown fee Fe/Cu(100) ultrathin films, 
as the temperature decreases from T c , even though those 
systems do not present SRT—. 



In the Saturated state, the stripe width increase with 
rj is directly related to the wall width decrease through 
the relation h ~ e s / 2l ^_ The wall width in turns is 
determined by the competition between exchange and 
anisotropy. Once the anisotropy is large enough that the 
wall width reaches the atomic limit, h growth stops. One 
may wonder whether a similar mechanism could be be- 
hind the stripe width variation with temperature. 

Besides its direct application to real systems, know- 
ing the ground state of this system for arbitrary values 
of the exchange coupling is of fundamental importance 
to have a correct interpretation of Monte Carlo simula- 
tion results. Being one of the more powerful tools to 
analyze these kind of systems at the present (specially 
at finite temperatures), it is basically limited by finite 
size restrictions, which implies relatively small values of 
S (the characteristic length h of the problem grows expo- 
nentially with 8 at low temperatures). 
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APPENDIX A: YAFET & GYORGY 
APPROXIMATION 

We briefly review in this appendix the derivation of 
the main results of YG approximation 1 ^. 

1. Energy per spin (Eq.©) 

The expression for the exchange and anisotropy ener- 
gies per spin from Eq.JT]) in terms of the one dimensional 
magnetization profile M(x) is straightforward: 

e cx + e an = S-^J2 M(x).M(x + 1) - J i MZ ( x )f 

X X 

(Al) 

The dipolar energy per spin can be expressed as e<np = 
e d.ip + e d%> w h ere e di P ls the self-interaction energy (i.e., 
the sum over x of the interaction energy between spins 
belonging to the same line at x) and e^™* is the interac- 
tion energy between all different pairs of lines. The self 
interaction term is given by^ 

e s dip = -2g + ^J2l MZ ( x )? ( A2 ) 

X 

with 
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oo 

* = E=s = c(3) = i. 
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(A3) 



where the sum in the above expression is taken over all 
values of (x, x') such that x ^ x' . The interaction energy 
between two lines located at x and x + n is given by 



where ((x) is the Riemann Zeta function. The interaction 
term can be expressed as 

^p^E^m ( A4 ) 

x^x' 



Ef nt (x, x + n) = M z (x) M z (x + n)f x {n) + M y (x) M v (x + n) - 3 / a (n)) 

I 



(A5) 



where 



where ci = /i(l) — 1 and c-i = /i(l) — 3/2(1). Finally 



2L ^ [n 2 + (v - y') 2 f 2 



m») = ^E 



(y y') 



'\2 



[n 2 + (y - y'ff 2 



(A6) 



(A7) 



In the limit L — > 00 the sums in Eqs. l|A6| and (| A7[) 
can be evaluated using a continuum aproximation^ giv- 
ing fi(n) ~ 1/n 2 and /i(n) — 3/2 (n) ~ 0. For n > 1 
the error in this approximation is smaller than 0.1%. 
For n = 1 they can be evaluated numerically giving 
/i(l) = 1.01243... and /i(l) - 3/ 2 (l) = 0.07276.... 
Then, Eq . (|A5[) can be written as 



int 1 ^ M*(x)M*(x') 
" dvp L t-i. \x-x'\ 2 



Xj^x' 

2ci J E M{x).M(x + 1) + C [M y {x)\ (A9) 



where C[M y (x)] is given by Eq.|f4]). Putting all these 
terms together we get Eq.(fS]). 



2. Variational equations for a striped 
magnetization profile with sinusoidal wall (SWP) in 
the continuum limit 



M z (x) M z {x + 1) + ci M(x).M(x + 1) 
(c 2 -c 1 )M y (x)M v (x + l). (A8) 



In the continuum limit h ^> 1 and is > 1 (5 > 1) 
Eq. lfl4|) can be written as 



1 1 k M 2 

M,S,rj =S-5-J2co S (Hx)-Hx + l)) + M 2 £ 6 ^/j TO (ft) _ _^ £ {£, (A10) 



m— 1,3, 



m— 1,3,... 



where k = 77— 3 <?+2(c2 — c\) and the functions D m (h) are 
given in Eq. (fT0|) . From Eqs. lfT5|) the Fourier cofficients b m 
in this limit are given byi£ 



b m = (-l)^- 1 )' 2 — l —cos ( ^) (All) 



ancU2 



E bl = (2-A) 



m— 1,3, . 



(A12) 



with A = w/h. The dipolar energy term in Eq. ljAlOj) can 
be approached byi^ 



£ b 2 n D m (h) = 

m— 1,3,. . . 



?L (2 _ A) _Z^ G(A ) 



(A13) 
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where we have used Eq. (|A12p . the cuadratic term in Taking cos(<p(x)) — cos(9) cos(irx/w) in the region inside 



Eq. lflOl) has been neglected, k = ir/h and 

G(A) = J2 mb "n( A ) 
m=l ,3, . . . 



(A14) 



Assuming a smooth profile inside the walls \<fi(x) — (p(x + 
1)| <C 1 the exchange term in Eq. (|Al0p can be ap- 
proached by 



a wall and replacing the summation in Eq. (|Al5P by an 
integral we have that 



S+^(l-ainO) (A16) 



2w ^-i I dx 

x=l v 



Replacing Eqs. (|Al2|l . I|A13|) and (|A16|) into Eq. (|AlO|l 
(A15) we get 



&swp 



2A 



k\l 



+ 7 1 



7T k 



G(A) cos 2 e 



(A17) 



r 



where 7 = 7r 2 /3 — k. Minimizing Eq. (|A17P respect to the 
variational parameters (9, A, k) we get 



5k 
A" 

5_k 2 

A 



(1 



— (1-sint 



= -G(A)cos 2 9 
, dG 

J + TTK — — 

dA 



(A18) 
(A19) 



Both G(A) and the expression between square brackets 
in Eq. (|A23P are monotonously decreasing functions of 13 
A. Since the maximum allowed value is A = 1 (which 
corresponds to w = h. i.e., pure sinusoidal profile), the 
minimum value of k = K, m in for which a domain solution 
exists corresponds to A = 1. Using thali^ 5 - G(l) — 1 and 
(dG/dA)A=i = — 1 we have 



Sk 2 
~2A 



2n 



l--)-nkG(A) 



sin a cos t 



(A20) 

Notice that The fully saturated state 9 = is never a 
solution of the above equations, except in the limit k — > 
(or A — > 0), which corresponds to S — > 00. On the other 
hand, the planar ferromagnetic state 9 = it/2 (e = 0), is 
always solution of the above equations. For 9 ^ n/2, the 
variational equations reduce to Eqs ((T7|l . 

Close to the SRT (i.e., the transition between a state 
with 9 = tt/2 and one with 9 ^ n/2) we can assume 
cos 6 = s <C 1 and therefore 1 — sin 9 = 1 — vl — s 2 ~ 
s 2 /2. Replacing into Eqs. (jA18[) - (|A20j) they become 



k 

Sk 2 
2A2 



7rA 



G(A) 

, dG 

7 + 7T k — — 

dA 



(A21) 
(A22) 



independent of s, while Eq. i|A20|) becomes identically 
zer o in th e limit s — > (SRT). Replacing Eq. (|A21|l into 
Eq.{A22| we find 



7T- 



k _ 1 G(A) 
~ ~ 3 + 



25 



G ,A, + 2Af 



(A23) 



Kmin -I |j M 



(A24) 



From Eq. (|A2ip this corresponds to fc = 7r/<5 or h = S. 
Replacing these values into Eq. (TT8|l we see that in the 
limit k — > K m in we have that sin 9 — > 1. Also from 
Eq . (j A17[) we see that in this limit the planar ferro- 
magnetic and the domain solutions become degenerated. 
Hence, this point corresponds to the SRT and the SRT 
line in the (n, 5) space is given by Eq. (fl9|) . 



APPENDIX B: EXACT ENERGY EVALUATION 
FOR A STRIPED MAGNETIZATION PROFILE 
WITH SINUSOIDAL WALL (SWP) IN A 
LATTICE 

It is easier to carry out this calculation by considering 
a profile whose wall starts at x = 0, i. e. 



M{x) 



M cos (2^) if 0<x<w 
—Mq if w < x < h 



(Bl) 



where h = 1,2,... and w = 1,2,..., ft,. The pro- 
file Eq. (fl5|) is related to the previous one by M z (x) = 
M z (x - (h - w)/2). The profile (|Bl]l can be expanded 
as 
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M{x) 



Mo 
2h 



2/1+1 

E 

m. — 1,3.. . 



(B2) 



21) 



K 



C m € 



ii^ra (h — w) jlh 



(B3) 



where lZ[z] stands for the real part of z. The coefficients 
so the coefficients b m of the expansion J7]) are given by c m are given by 



1 2h w—1 h—l—w 

Cm = Jf Yl M z (x)e~ iS ir x = 2 V cos (—) e~ l ^ x - 2 V <-'^> r+ "' ' 



'0 



:r=0 



x=0 



(B4) 



The summations in Eq. (|B4|l can be carry out explicitly 
obtaining 



APPENDIX C: ZERO TEMPERATURE MONTE 
CARLO TECHNIQUE FOR STRIPED DOMAIN 
PATTERNS 



= f + + f ~ + f° 



(B5) 



where 



x=0 



1-e 
ID 



if, 



^1 



if e ITrQ± =1 



with a± = (±i - f ) and 



h—l — w 

2 g-jxW 



= 2 



^iirmw j h _|_ ^ 
g—iirmfh 1 



(B7) 



The dipolar energy can then be evaluated from Eqs.Q 
and (fTUl) . The anisotropy energy can be easily calculated 
and gives 



-kM 2 



if w = 1 



-kM 2 (1-|) otherwise 



(B8) 



The exchange and correction terms in Eq. (fl4|) can be 

expressed as 



In order to check the different striped profiles used to 
minimize the energy of the system, we implemented a 
simulated annealing protocol, based on Metropolis dy- 
namics. The temperature was decreased down to zero at 
a constant rate T(t) = Tq — rt, where the time is mea- 
sured in Monte Carlo Steps. All the simulations are made 
starting from a planar ferromagnetic state with T = 1 
and r = 10~ 4 . Every simulation is repeated 100 times us- 
ing a different sequence of random numbers, to check for 
possible trapping in local minima. Since we are consider- 
ing only periodic straight domains with Bloch walls, the 
problem is basically one dimensional and we can restrict 
the search to a one dimensional pattern over the x direc- 
tion fixing M x {x) = and imposing periodic boundary 
conditions (PBC) in the y direction. We also use PBC in 
the x direction. In other words, we simulated a a lattice 
with L x x L y with L y = 1 and PBC, which are imple- 
mented by means of the Ewald sums technique. For every 
set of values of (<5, 77) we check the results for different val- 
ues of L x in order to avoid artificial frustration. We also 
performed some comparisons with MC results in a square 
Lx — Ly lattice; the results were indistinguishable. This 
ansatz allows us to obtain MC results for values of 5 up 
to 5 = 8 (for which the maximum equilibrium value is 
ft = 48). 



, w-l APPENDIX D: FOURIER COEFFICIENTS FOR 

e exc = -5'~ X M(x).M(x + 1) - S'(l - A) M 2 (B9) THE HYPERBOLIC TANGENT WALL PROFILE 



h x =o ' (HWP) 

The function tanh(z) is very well approximated by 
&(!- x) if 0<x<± 



C [M v {x)] = {C2 h Cl> X M y {x)M y {x + 1) 

+ 2(c 2 - ci)(l - A) (1 - M 2 ) (BIO) tanh(2:) ~ 1 C 1 - e-~ 2 *) 2 (l + e^) if \ < x 

(Dl) 

The summations in the above equations involve a finite Then the Fourier coefficients b m for the profile Eq. (|20|) 
number of terms that can be computed explicitly. can be expressed as b m — b] n + b 2 m) where 
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where h r = h/l w . Both integrals can be solved ana- 
x 3 lytically, leading to rather long expressions that can be 

b i = ± I x (i _ sm ( m7TX j ^ x (T)2) handled with symbolic manipulation programs. 
h r Jq 3 \ h r J 



and 



t •»=H^ (l ^ )2(l+e " I)sm ( ! ^)' ,I 

(D3) 
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